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ABSTRACT 

We consider the general expressions for the time delay of photons of different 
energies in the framework of multi-dimensional cosmological models. In models 
with compactified extra-dimensions (Kaluza-Klein type models), the main source 
of the photon time delay is the time variation of the electromagnetic coupling, 
due to dimensional reduction, which induces an energy-dependence of the speed 
of light. A similar relation between the fine structure constant and the multi- 
dimensional gauge couplings also appears in models with large (non-compactified) 
extra-dimensions. For photons of energies around 1 TeV, propagating on cosmo- 
logical distances in an expanding Universe, the time delay could range from a 
few seconds in the case of Kaluza-Klein models to a few days for models with 
large extra-dimensions. As a consequence of the multi-dimensional effects, the 
intrinsic time profiles at the emitter rest frame differ from the detected time pro- 
files. The formalism developed in the present paper allows the transformation of 
the predicted light curves of various energy ranges of the emitter into the frame 
of the observer, for comparison with observations. Therefore the study of energy 
and redshift dependence of the time delay of photons, emitted by astrophysi- 
cal sources at cosmological distances, could discriminate between the different 
multi-dimensional models and/or quantum gravity effects. 

Subject headings: cosmology-extra-dimensions: gamma rays: bursts-radiation 
mechanisms: photon delay 



1. Introduction 



One of the most challenging issues of modern physics is the existence of the extra- 
dimensions, idea proposed originally by Kaluza (1921) and developed by Klein (1926). Multi- 
dimensional geometries are the natural framework for the modern string/M theories (Witten 
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1996) or brane models (Horava and Witten 1996). String models also provide a natural and 
self-consistent explanation for the possible variation of the fundamental constants, as initially 
suggested by Dirac (1937, 1938, 1979). Hence the problem of the extra-dimensions of the 
space-time continuum is closely related to the problem of the variations of the fundamental 
constants, like, for example the fine structure constant or the speed of light (for a recent 
review of the experimental and theoretical studies and the present status of these fields see 
Uzan (2003) and Magueijo (2003). Most of the theories with extra-dimensions contain a 
built-in mechanism, which allows the variation of the fundamental constants. Within the 
multi-dimensional approach the physical interactions are described by a theory formulated 
in 4 + D dimensions, and the conventional four-dimensional theory appears as a result 
of a process of dimensional reduction. Couplings in four dimensions are determined by 
a set of constants of the multidimensional theory and the size A of the space of extra- 
dimensions. The multi-dimensional constants are assumed to be genuinely fundamental and, 
consequently, they do not vary with time. On the other hand it is natural to assume that 
in an astrophysical or cosmological context A varies with time, similarly to the scale factor 
a of our four-dimensional Universe. But this leads to the time variation in four dimensions 
of the fundamental constants, like fine structure constant a or the gravitational coupling G. 
Moreover, since their time dependence is given by the same factor A, the time variations of 
a and G could be correlated (Langacker et al. 2002). 

The search for a unification of quantum mechanics and gravity is likely to require a dras- 
tic modification of the present day deterministic representation of the space-time properties. 
There is at present no complete mathematical model for quantum gravity, and no one of the 
many different models proposed so far can give a satisfactory description of the physics on 
characteristic scales near the Planck length lp. However, in several of the approaches trying 
to find a theory of quantum gravity the vacuum can acquire non-trivial optical properties, 
because of the gravitational recoil effects induced by the motion of the energetic particles. 
The recoil effects may induce a non-trivial refractive index, with photons at different ener- 
gies travelling at different velocities (Ellis et al. 2000a). Photon polarization in a quantum 
space-time may also induce birefringence (Gambini and Pullin 1999), while stochastic effects 
in the vacuum could give rise to an energy dependent diffusive spread in the velocities of 
different photons (Ellis et al. 2000b,c). 

Therefore a large class of physical models, incorporating quantum gravitation and/or 
multi-dimensional field theories predict that the propagation of particle in vacuum is mod- 
ified due to the supplementary effects induced by the modification of the standard general 
relativity. In particular, the possible violation of the Lorentz invariance or the existence of 
extra-dimensions can be investigated by studying the propagation of high energy photons 
emitted by distant astrophysical sources (Amelino-Camelia et al. 1998). 
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The highest energy extra-galactic 7-ray sources in the known universe are the active 
galaxies called blazars: objects that emit jets of relativistic plasma aimed directly at us. 
Objects known as high frequency BL Lac objects (HBLs) are expected to emit photons in 
the multi-TeV energy range. Only the nearest ones are expected to be observable in TeV 
energies, the others being hidden by intergalactic absorption (Stecker et al. 1996). 

Extragalactic photons with the highest energy yet observed originated in a powerful flare 
coming from the giant elliptical active galaxy known as Markarian 501 in 1997 (Aharonanian 
et al. 1997) and from Markarian 421 in 2001 (Krennrich et al. 2001). The high energy flux of 
these emissions permitted detailed spectra to be extracted. Since as many as 25, 000 photons 
were detected, the spectra were derived with a high statistical accuracy. The observation 
of the TeV photons from the Markarian 501 and 421 allows to impose some constraints 
on the quantum gravity scale and on the breaking of the Lorentz invariance (Biller et al. 
1999; Stecker 2003). Quantum gravity phenomena are a result of the quantum fluctuations 
on the Planck scale Mp = ^hc/G. In models involving large extra-dimensions, like, for 
example in the brane-world models (Randall and Sundrum 1999a, b) the energy scale at 
which gravity becomes strong can be much smaller than Mp, with the quantum gravity 
scale Mq G approaching the TeV scale. 

The data from the TeV gamma -ray flare associated with the active galaxy Markarian 421 
have been used to place bounds on the possible energy dependence of the speed of light in the 
context of an effective quantum gravitational energy scale in (Biller et al. 1999). The limits 
derived indicate this energy scale to be higher than 6 x 10 16 GeV for the approach to quantum 
gravity in the context of D-brane string theory. By assuming a modified dispersion relation 
between the energy E 7 and the momentum p 7 of the photon of the form = j£ f — p^/Mqg, 
and a similar relation for the electron, one can obtain the constraint M QG > E^/8ml for the 
quantum gravity scale (Stecker 2003). Since pair production occurs for energies of at least 
20 TeV, it follows that M QG > 0.3M P . The results also indicate an absence of evidence for 
the violation of the Lorentz invariance, as proposed by some quantum gravity and multi- 
dimensional models. 

Strong constraints on Lorentz violating microscopic structures of space-time, like, for 
example, discreteness, non-commutativity, or extra dimensions can be obtained from the 
observation of 100 MeV synchrotron radiation from the Crab nebula (Jacobson et al. 2003). 
The Crab synchrotron emission has been observed to extend at least up to energies of 100 
MeV, just before the Compton hump begins to contribute to the spectrum (Atoyan and 
Aharonian 1996; Hillas et al. 1998). The magnetic field in the emission region has been 
estimated to a value between 0.15 — 0.6 mG (Hillas et al. 1998). To produce the observed 
radiation of 100 MeV in this field requires a relativistic 7-factor of the order of 7 = 3 x 10 9 , 
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corresponding to an electron energy of the order of 1500 TeV, with an electron velocity 
differing from c by less then 10~ 19 c. Then the observation of 100 MeV synchrotron radiation 
from the Crab nebula gives the constraint Eqg > 10 26 GeV (Jacobson et al. 2003). Hence 
this observation rules out this type of Lorentz violation, providing an important constraint 
on theories of quantum gravity and imposes a stringent constraint on any modification of 
the dispersion relations of the electron that might be induced by quantum gravity. 

The spectrum of the synchrotron radiation from the coupling of an electrically-charged 
particle to an external magnetic field, in the presence of quantum-gravity effects of the 
general form (E/Mqg)" has been derived in Ellis et al. (2003a). The synchrotron constraint 
from the Crab Nebula practically excludes a ~ 1.74 for M QG ~ M P ~ 1.2 x 10 19 GeV. 
The model suggests a linear modification of the dispersion relation for the photon, but not 
for the electron, and hence is compatible with known constraints from the Crab Nebula. 
New constraints on possible Lorentz symmetry violation of the order of E/Mp for electrons 
and photons in the framework of effective field theory, by using the absence of vacuum 
birefringence in the recently observed polarization of MeV emission from a gamma ray burst, 
and the absence of vacuum Cerenkov radiation from the synchrotron electrons in the Crab 
nebula, have been derived in Jacobson et al. (2003). These constraints allow to improve the 
previous bounds by eleven and four orders of magnitude, respectively. 

The possibility of the use of the high energy radiation from gamma-ray pulsars to place 
limits on quantum gravity effects has been suggested by Kaaret (1999). The emission from 
the Crab pulsar at energies above 2 GeV trails that at 70 — 100 MeV by no more than 
0.35 ms (95% confidence). This effect places a lower bound on the energy scale of quantum 
gravitational effects on the speed of light of 1.8 x 10 15 GeV. In the near future this bound 
might be improved by two orders of magnitude by observation of pulsations from the Crab 
at higher energies, of the order of 50 — 100 GeV. 

The confirmation that at least some gamma-ray bursts (GRBs) are indeed at cosmo- 
logical distances raises the possibility that observations of these could provide interesting 
constraints on the fundamental laws of physics (for recent reviews on GRBs see (Zhang 
and Meszaros 2003; Cheng and Lu 2001)). The fine-scale time structure and hard spectra of 
GRB emissions are very sensitive to the possible dispersion of electromagnetic waves in vacuo, 
with velocity differences Am ~ E / Eqg, as suggested in some approaches to quantum gravity. 
GRB measurements might be sensitive to a dispersion scale Eqg comparable to the Planck 
energy scale Ep ~ 10 19 GeV, sufficient to test some of these theories (Amelino-Camelia et 
al. 1998). Hence the study of short-duration photon bursts propagating over cosmological 
distances is the most promising way to probe the quantum gravitational and/or the effects 
related to the existence of extra-dimensions. The modification of the group velocity of the 
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photons by the quantum effects would affect the simultaneity of the arrival times of photons 
with different energies. Thus, given a distant, transient source of photons, one could mea- 
sure the differences in the arrival times of sharp transitions in the signals in different energy 
bands. A key issue in such a probe is to distinguish the effects of the quantum-gravity /multi- 
dimensional medium from any intrinsic delay in the emission of particles of different energies 
by the source. The quantum-gravity effects should increase with the redshift of the source, 
whereas source effects would be independent of the redshift in the absence of any cosmologi- 
cal evolution effects. Therefore it is preferable to use transient sources with a known spread 
in the redshift z. The best way to probe the time lags that might arise from quantum gravity 
effects is to use GRBs with known redshifts, which range up to z ~ 5. 

Data on GRBs may be used to set limits on variations in the velocity of light. This has 
been illustrated, by using BATSE and OSSE observations of the GRBs that have recently 
been identified optically, and for which precise redshifts are available, in Ellis et al. (2000d). 
A regression analysis can be performed to look for an energy-dependent effect that should 
correlate with redshift. The analysis of GRBs data yield a limit Mqg ~ 10 15 GeV for the 
quantum gravity scale. The study of the the times of flight of radiation from gamma-ray 
bursts with known redshifts has been considerably improved by using a wavelet shrinkage 
procedure for noise removal and a wavelet 'zoom' technique to define with high accuracy 
the timings of sharp transitions in GRB light curves (Ellis et al. 2003b). This procedure 
optimizes the sensitivity of experimental probes of any energy dependence of the velocity of 
light. These wavelet techniques have been applied to 64 ms and TTE data from BATSE, and 
also to OSSE data. A search for time lags between sharp transients in GRBs light curves in 
different energy bands yields the lower limit Mqg > 6.9 x 10 15 GeV on the quantum-gravity 
scale in any model with a linear dependence of the velocity of light, c ~ E /Mqg- 

High energy GeV emissions from GRBs have already been detected (Sommer et al. 1994; 
Hurley 1994). There is also tentative evidence for TeV emissions (Atkins et al. 2000; Poirier 
et al. 2003). The production of TeV photons is also predicted by most of the GRB theories. 
The emission mechanisms for TeV photons include electron inverse Compton emission and 
synchrotron emission from the protons accelerated by GRB shocks. The shocks could be 
internal shocks, external forward shocks or external reverse shocks of GRBs. Such very high 
energy photons at cosmological distances may largely be absorbed by interacting with the 
cosmic infrared background radiation (Madau and Phinney 1996). 

There are many astrophysical mechanisms that could produce a delay in the arrival 
time of high energy photons or cosmic rays. The electron inverse scattering of the created 
electron-positron pairs off the cosmic microwave background photons will produce delayed 
MeV-GeV emission. There are two likely mechanisms causing the time delay. One is the 
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angular spreading effect of the secondary pairs, that is, the scattered microwave photons 
deviate from the direction of the original TeV photons by an angle ~ I/7, where 7 is the 
Lorentz factor of the electron-positron pairs (Cheng and Cheng 1996; Dai and Lu 2002; Wang 
et al. 2004). Another mechanism is related to the deflection of the direction of propagation 
of the pairs in the intergalactic magnetic field, if this field is sufficiently strong (Plaga 1995). 
It is important to note that all these delay mechanisms predict that the low energy photon 
come late. But in the models with extra-dimensions the high energy photon comes later. 

It is the purpose of the present paper to point out some other possible sources of time de- 
lay for high energy photons from GRB emissions. Namely, we shall consider the effects of the 
existence of the extra-dimensions on the propagation of high energy photons in two distinct 
physical scenarios: Kaluza-Klein theories with compactified extra- dimensions and models 
with large extra-dimensions. In Kaluza-Klein theories the time variation of the scale factor 
of the space of extra dimensions leads to a time variation of both the gravitational constant 
G and of the fine structure constant a. The variation of a is photon-energy dependent, and 
via the time delay allows the estimation of the size of the (compactified) extra-dimension. 
In the case of the large (non-compactified) extra-dimensions, there is no multi-dimensional 
physical mechanism to induce a variation of a. Recently, using high resolution spectroscopy 
of QSO absorption spectra, a time variation of a has been reported (Webb et al. 1999; Mur- 
phy et al. 2001, 2003). The detected rate of change of a is of the order Aa/a ~ — 10~ 5 at a 
redshift z ~ 1.5. 

The present paper is organized as follows. In Section II we discuss the most important 
physical processes that could lead to the variation of the electromagnetic coupling, in the 
different versions of the mult i- dimensional models. The basic equations for the time delay 
of photons in Kaluza-Klein and Randall-Sundrum type models are obtained in Section III. 
In Section IV we discuss and conclude our results. 



2. Variation of the electromagnetic coupling in models with extra dimensions 

The starting point in the multi-dimensional approach of describing the fundamental 
interactions by Kaluza and Klein is the consideration of the pure Einstein gravity in a 
multi-dimensional space-time M 4 x K D , described by the multi-dimensional metric tensor 
cjmn- Here M 4 is the four-dimensional space-time and K D is the compact manifold of 
extra-dimensions. Generally the reduced theory contains the Einstein gravity and the Yang- 
Mills fields with the gauge group determined by the isometry group of the space of extra- 
dimensions. The action of the multi-dimensional Kaluza-Klein theory is the action for the 
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pure Einstein theory on M 4 x K D , with the action given by (Overduin and Wesson 1997) 



I [ d 4+D x^R( 4+D \ (1) 

t(4+D) J 



1 (i;; t 

where g = det (cimn)i R^ +L)S> is the scalar curvature in M 4 x K D and G( 4+ _d) is the multi- 
dimensional gravitational constant, which is assumed to be a true constant and does not 
depend on time. According to the procedure of dimensional reduction, the (//, v — 0, 1, 2, 3) 
components of the metric tensor are identified as the four-dimensional metric tensor, while 
certain combinations of the rest of the components are identified as gauge field multiplets 
and scalar fields m , m — 1,2.... After the mode expansion of these fields, with the coefficients 
of the expansion (interpreted as four-dimensional fields) depending only on x M , /i — 0, 1, 2, 3 
we obtain the four-dimensional action 



d 4 x 



\ p(4) i V —Tr (F^ F^) 



(2) 



where G(t) = G<±(t) is the four-dimensional gravitational constant. In obtaining Eq. (2) we 
have considered only zero modes of the mode expansion. The parameters gi (t) are the gauge 
couplings, and the index i labels the simple subgroups of the gauge group. The general 
reduction of the initial Kaluza-Klein action S also gives terms including non-zero modes of 
the gravitational, gauge and scalar fields (Overduin and Wesson 1997). The scalar fields 
give highly non-linear interaction terms and are coupled non-minimally to the gravitational 
and gauge fields. In the following we neglect, for simplicity, the contribution of these scalar 
fields. 

Identifying the gravitational and gauge couplings from the action 5o for the zero modes 
one obtains the following expressions for G(t) and gfit) in terms of G(4+d) and the radius 
$ (t) of the space of the extra-dimensions (Loren-Aguilar et al. 2003): 

a (t) - G(4+D) n 2 (t) - k G(4+D) m 

where Vd (t) ~ $ D (t) is the volume of the space of the extra-dimensions and hi are coeffi- 
cients which depend on the isometry group of K D . The fine structure constant a (t) is given 
by a linear combinations of gf{t), with the specific relation depending on the gauge group 
and the scheme of spontaneous symmetry breaking. Generally, 

where k\ is a constant. Since Vd (t) /Vd (t) = d (^/^j, for the time variation of the fine 
structure constant we obtain 

a $ 

- = -P + 2) 5 - (5) 
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The variation of the fine structure constant in a more general model, in which a Yang- 
Mills type field is also included in the (4 + D)-dimensional space-time was also considered 
in Loren-Aguilar et al. (2003). The action is 



/ 



d i+D x\f— 



9 



1 



16nG 



(4+D) 



% 2 



-Tr 



(4+D) 



(f mn f mn ) 



(6) 



with g 2 i+D \ the mult i- dimensional coupling, supposed to be constant in time. In this case the 
dimensionally reduced theory includes the Einstein gravity and the four-dimensional gauge 
fields plus scalar fields with a quartic potential. The four dimensional gravitational and fine 
structure constants are given by 



G(t) 



G 



(4+D) 



a (t) = h 



G 



(4+D) 



" z Vb(t)' 

where k 2 is a constant. The time variation of the fine-structure constant is 



(7) 



a $ 
- = -D-. 

a $ 



(8) 



We consider now models with large extra-dimensions (brane world models), as initially 
considered in Randall and Sundrum (1999a). In the 5D space-time the brane-world is located 
at Y (X 1 ) = 0, where X 1 , / = 0, 1, 2, 3, 4, are 5-dimensional coordinates. The effective action 
in five dimensions is (Maeda and Wands 2000; Chen et al. 2002) 

S = J d 5 X^F 5 - A 5 ) + J y Q d'x^-g (±K± - A + L matter ) , (9) 

where = 8nG 5 the 5-dimensional gravitational coupling constant, A 5 is the cosmo logical 
constant in the bulk. x M , /i — 0, 1,2,3, are the induced 4-dimensional brane world coordi- 
nates. -R5 is the 5D intrinsic curvature in the bulk and K ± is the extrinsic curvature on 
either side of the brane. 

Assuming a metric of the form ds 2 = (rijrij + gi^dx 1 ' dx 3 ' , with riidx 1 = d\ the unit 
normal to the x — constant hypersurfaces and gu the induced metric on x — constant 
hypersurfaces, the effective four-dimensional gravitational equations on the brane take the 
form (Shiromizu et al. 2000; Sasaki et al. 2000; Chen et al. 2001a,b): 

G>„ = -A^ + kjT^ + k^S^ - E^, (10) 

where 

S, v = -TT, V - -T, a T ua + -g, u (3T^T a/3 - T 2 ) , (11) 
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and A = fc|(A 5 + /c 2 A 2 /6)/2, k\ = kjX/6. Ejj = CiA.jBn A n B . Ciajb is the 5-dimensional 
Weyl tensor in the bulk and A is the vacuum energy on the brane. T^ v is the matter energy- 
momentum tensor on the brane and T = T M ^ is the trace of the energy-momentum tensor. 

The reduction formula expressing the four-dimensional Planck mass Mpi in terms of 
the fundamental (five-dimensional) mass scale M = ( 167rG(5) j ~ k has been derived in 
Boos et al. (2002). The result is 

Ml, = ^ (e 2 ^ - 1) . (12) 

Since in the models with large extra-dimensions all matter fields are localized on the 
brane and do not depend on the radius of the extra-dimension, there is no simple mechanism 
for the variation of the fine structure constant. A possible form of the variation of a has 
been obtained in Loren-Aguilar et al. (2003), giving 

a{t) = h wr* = ~^ (13) 

where k% is a constant. 

Finally, we shall briefly consider the problem of the quantum corrections to the expres- 
sions for the fine structure constant presented above. This question has been discussed in 
Langacker et al. (2002) and Loren-Aguilar et al. (2003). To relate the value of a obtained 
at the scale M$ = (t) to its value at some low energy scale \x by taking into account 
quantum corrections one must use the renormalization group formulas for running couplings, 
which gives for the time variation of a 

^4 = - « (M) 4 (1 + M*) , (14) 

where A is a constant of order one. 

The second term in Eq. (14) is of order O (a). Hence it is a subdominant term, thus 
showing that quantum effects does not significantly affect the time variation of the fine 
structure constant. 



3. Time delay of photons in multi-dimensional expanding Universes 



We consider first the propagation of gamma-rays from GRBs in the Kaluza-Klein type 
models. For simplicity we restrict our discussion to the five-dimensional case. Hence we 
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assume a flat Friedmann-Robertson- Walker type background metric of the form 

ds 2 = c 2 dt 2 - a 2 (t) [dr 2 + r 2 {d6 2 + sin 2 6d V 2 )] + e$ 2 (t) dv 2 , (15) 

where a is the scale factor of the Universe, e — ±1 and $ is the scale factor of the fifth 
dimension, denoted by v. We also assume that the time variation of the fine-structure 
constant is entirely due to the change in the speed of light c. Therefore we neglect any 
possible time variation of the electric charge or Planck's constant. Then the time variation 
of the speed of light can be related to the size of the fifth dimension by means of the general 
equation 

Ac $ 

-^ = fe ( 16 ) 
Ac $ v ' 

where f3 — 1 in the case of the Einstein- Yang-Mills model and (3 = 3 for the case of the pure 

Einstein gravity in five dimensions. Eq. (16) can be integrated to give 

c = c (l+^), (17) 

where c is an arbitrary integration constant. In order to find a simple and directly testable 
relation between the radius of the extra-dimension and the energy of the photon, we shall 
assume, following the initial proposal in Ma (1990a,b) that the mass of a body (and the 
associated energy) corresponds to the length of a "line segment" of the fifth subspace. Such 
a relation embodies the spirit of the Mach's principle in the sense that the inertial mass 
depends on the distribution of the matter in the Universe. In a more general formulation, we 
shall assume that the variables parameters c, G and the photon energy E = hv are related 
to the metric tensor component of the fifth dimension by means of the equation (Mak and 
Harko 1999) 

G (t) E 1 



c 4 7 



y/\gu\dv = — §dv, 
7 Jv° 



where 7 is an arbitrary constant. 



If $ is independent of v, as is the case in models with compactified extra-dimensions, 
Eq. (18) gives 

jGE E_ 

~ c A {y — v°) ~ E K ' ( } 

where we denoted Ek = c^Av/'yG, with Av = v—v° describing the variation of the size of the 
fifth dimension between the moments of the emission and detection of a photon. Therefore 
the energy-dependence of the speed of light of the photon is given by 



c 



(20) 



- 11 - 



We consider two photons emitted during a gamma-ray burst with present day energies 
Ei and E 2 . At earlier epochs, their energies would have been blue-shifted by a factor 1 + z. 
Then it follows that the difference in the velocities of the two photons is given by 

Ac = Co AE£ ±£ ) /3=li 

£j K 

and 

Ac =t /<*-^< 1 + «>' ,/» = 8, (22) 

respectively, where we denoted = E x — E 2 and / i? 2 ) = Ef — E\. A linear energy 
dependence of the difference of the photon velocities has also been considered in (Ellis et 
al. 2003b) as a result of the dispersion-relation analysis of the Maxwell equations in the 
non-trivial background metric perturbed by the recoil of a massive space-time defect during 
the scattering of a low energy photon or neutrino. 

For light propagating from cosmological distances the differential relation between time 
and redshift is (Ellis et al. 2003b) 

dt = ~ H °\l + z)h(zY (23) 

where 



g(z) = ^n A + n M (l + z) 3 , (24) 

Ho = 72 km s" 1 Mpc -1 (Freedman et al. 2001), Qm ~ 0.3 and Q\ 0.7 are the mass density 
parameter and the dark energy parameter, respectively (Peebles and Ratra 2003). A particle 
with a velocity c travels an elementary distance cdt = —H 1 cdz/ (1 + z) g(z), with the differ- 
ence AL in the distances covered by the two particles given by AL = Hq 1 J q z Acdzj (1 + z) g(z) 
(Ellis et al. 2003b). By taking into account the expression for Ac we obtain the following 
equations describing the time delay of two photons: 

%P = h (25) 



At 


= H ^ E 

IliK 


lo 


Ho 1 


\f(Ei,E 2 y 

E K 


7 

J 



'1 + z) 2 dz 



9(z) 



,(3 = 3. (26) 



In the case of isotropic homogeneous cosmological models with large extra-dimensions 
there is a non-zero contribution from the 5-dimensional Weyl tensor from the bulk, expressed 
by a scalar term U, called dark radiation (Chen et al. 2002; Harko and Mak 2003; Chen et 
al. 2003). The "dark radiation" term is a pure bulk (five dimensional) effect, therefore we 
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cannot determine its expression without solving the complete system of field equations in 
5 dimensions. In the case of a Friedmann-Roberston- Walker type cosmological model the 
expression of the dark radiation is (Chen et al. 2001a; Dabrowski et al. 2002) 

U=% (27) 
a 4 

with Uq an arbitrary constant of integration. Since the fifth dimension is large, the scale 
factor $ can also be a function of v. Hence an explicit knowledge of the v-dependence of $ is 
needed in order to derive the speed of light- photon energy dependence. However, taking into 
account Eq. (13), which shows a linear dependence of a on the scale of the fifth dimension, 
we can assume that the time delay between two different energy photons emitted during a 
gamma ray burst is given by 

_ t AE f z dz 



At = H ° E7j WY (28) 

where Ep is the energy scale associated with the large extra-dimensions and 

h(z) = ^Jn A + n M (l + zf + flu (1 + z)\ (29) 
with VLjj the dark radiation parameter. 



4. Discussions and final remarks 

In order to calculate the delay in the gamma ray photons arrival time in the Kaluza- 
Klein type models we need to estimate first the Kaluza-Klein energy scale Ek = c A Av/^fG for 
which the effects from the extra-dimensions become important. Assuming that the size and 
the variation of the extra-dimension Av is of the same order as the Planck length, Av = lp = 
1.6 x 10~ 33 cm, the Kaluza-Klein energy scale is equal, for 7 = 1, to Ek = 1-2 x 10 19 GeV. Of 
course, a large value of 7 can decrease the Kaluza-Klein energy scale. The variation of the 
difference in the photon arrival time, for different values of the photon energy is presented, 
for the Einstein- Yang-Mills model corresponding to f3 — 1 and with 7 = 1, in Fig. I. 

Due to the power three energy dependence of the time delay in the pure Einstein gravity 
model (P = 3) the value of At is extremely small for 7 = 1, corresponding to the E K = 
1.2 x 10 19 GeV energy scale. In this case At ps 10~ 30 s even for photon energies of the order 
of 1 TeV. In order to obtain some observable effects a very large value of 7, of the order of 
7 = 10 10 is required. The variation of the photon time delay in the case of the Kaluza-Klein 
compactified model with pure Einstein gravity is represented in Fig. 2. 

The analysis of the BATSE and OSSE data has imposed a lower limit E QG > 6.9 x 10 15 
GeV on the quantum gravity scale in the linear model (Ellis et al. 2003b), which is much 
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smaller than the Kaluza-Klein energy scale we have considered. If the fundamental energy 
scale is of the order of Eq G) then the time delay between TeV and KeV/MeV photons could 
have larger values than those considered in the present approach. 

To generate the correct hierarchy between the Planck scale and the TeV scale in models 
with large extra-dimensions, the product /c$ must be of the order of k& « 11 — 12 (Boos 
et al. 2002). There are no simple mechanisms to describe the temporal evolution of the fine 
structure constant in this type of models, a variation of a requiring the consideration of bulk 
gauge and, perhaps, fermionic fields. However, based on the analogy with the Kaluza-Klein 
case one can assume a linear dependence of the time delay on AE, with the characteristic 
energy scale Ep a parameter to be determined from observations. The variation of the time 
delay for two different energy photons is represented, for E F = 7 x 10 15 GeV, in Fig. 3. 

In deriving the equations for the time delay due to the presence of the effects of the 
extra- dimensions we have used a crucial assumption, namely, we have considered that there 
is an intrinsic synchronization of pulsing emission of photons at different energy ranges. This 
assumption allows the consistent determination, by comparison with observations, of the pa- 
rameters determining the effects of the extra-dimensions on the photon propagation. How- 
ever, if the emission of the photons of different energies is not synchronized, the simple pro- 
cessing of data for the delay of the photons is not enough to constrain the multi-dimensional 
effects. There are many proposed explanations, which could produce a delayed time scale 
at the emitter between, for example, GeV and keV/MeV photons, as observed in the case 
of GRB 940217 (Hurley 1994). These explanations include interaction of TeV photons with 
cosmic infrared background photons (Plaga 1995), interactions of ultrarelativistic protons 
with a dense cloud (Katz 1994) or inverse Compton scattering in early forward and reverse 
shocks (Meszaros and Rees 1994). In these cases, in order to consistently determine the time 
delay between photons at different energies, due to multi-dimensional propagation effects, 
the knowledge of the initial time profiles of MeV, GeV and TeV photons is also required. 
The description of the non-synchronized emission time profiles is also model-dependent, and 
thus in this case it is difficult to clearly distinguish between emission and propagation effects. 
However, despite all these possibilities of producing a time delay of different energy photons 
at the emitter, in our model we make two very clear predictions, namely, (1) that the high 
energy (TeV) photons have the higher delay and (2) that the delay time scales are correlated 
with the redshift z of the emitting source, because the delay results from propagation effects. 
In the future, if enough observational data will be available, the study of the energy and red- 
shift dependence of the delay, At = At [E 1 , z) could lead to the possibility of discriminating 
between emission and propagation effects. 

We also suggest that if there is an intrinsic synchronization of the photon emission at 
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different energies, then due to the effects of the extra-dimensions the detected time profiles 
between the KeV/MeV/GeV and the TeV bursts should be very different. We want to 
emphasize that it is better to measure the time profile difference between KeV/MeV and 
TeV photons, instead of the difference between KeV/MeV and GeV photons, in order to 
avoid the contamination or other effects. For an emitter with an initial Gaussian time 
profile, e - (*/ r ) , where r is the duration of the burst, the pulse shape of the 1 — 10 TeV 
photons is shown in Fig. 4. We have considered that the GRB occurred at a redshift of 
z — 3, and considered the initial duration of the burst r = 1 s. For the time delay of the 
photons we have adopted the linear model, with the fundamental energy scale of the same 
order of magnitude as the Planck scale, E K = 1.2 x 10 19 GeV. In this case the time delays 
of the TeV photons with respect to the KeV/MeV/GeV photons are Ati = 52.9 s for the 1 
TeV photon and At 2 = 529.3 s for the 10 TeV photon, respectively. 

Realistically, in order to determine the exact amount of delay time, we have to measure 
the time profiles at the detector of both the KeV/MeV and of the TeV bursts, respectively. If 
the time delay between the KeV/MeV bursts and TeV bursts is longer than several seconds, 
the ground based TeV telescopes are capable to catch the burst. A good coordinate effort 
between the SWIFT satellite and the ground based TeV telescopes can easily make this 
measurement possible (Weekes 2003). 

Observations of the time delay in gamma ray bursts have been proposed up to now 
mainly as tests of quantum gravity effects. In a pioneering work Amelino-Camelia et al. 
(1998) suggested that quantum-gravitational effects could induce a deformed dispersion re- 
lation for photons of the form c 2 fP' = E 2 [1 + f {E/Eqg)\- By representing / in a form 
of a power series, the energy dependent velocity of the speed of light can be represented as 
v ~ c (1 ± E/Eq G ). This type of velocity dispersion results from a picture of the vacuum as 
a quantum gravitational 'medium', which responds differently to the propagation of parti- 
cles of different energies. This is analogous to propagation through a conventional medium, 
such as an electromagnetic plasma. The gravitational 'medium' is generally believed to 
contain microscopic quantum fluctuations, which may occur on scale sizes of order of the 
Planck length. In this approach the vacuum is viewed as a non-trivial medium containing 
'foamy' quantum-gravity fluctuations. The nature of this foamy vacuum may be visualized, 
for example, by imagining processes that include the pair creation of virtual black holes. 
The light propagation in the semiclassical space-time that emerges in canonical quantum 
gravity in the loop representation was considered in Gambini and Pullin (1999). In such a 
picture space-time exhibits a polymer-like structure at micro-scales, and departures from the 
perfect non-dispersiveness of an ordinary vacuum naturally occur. Maxwell equations are 
modified due to the quantum gravity, and non-vanishing corrections to the electromagnetic 
field equations appear that depend on the helicity of the propagating waves. These effects 
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could lead to constraints on the discrete nature of the quantum space-time from the study 
of gamma-ray bursts at different energies. 

In the present paper we have considered a different class of effects, which could generate 
a time delay At of the high energy photons, emitted during the gamma ray bursts, namely 
the possibility that the extra-dimensions of the space-time may modify the speed of light 
in vacuum. This is mainly due to the dependence of the fine structure constant a on the 
extra-dimensions. The effects related to extra-dimensions in the variation of a are much 
stronger than the quantum gravity effects. In models with extra-dimensions, the speed of 
light-photon energy dependence can be obtained exactly, by considering the background 
gravitational field and the scalar and Yang-Mills type gauge fields in extra-dimensions. In 
quantum-gravity models this dependence is modelled more or less phenomenologically. On 
the other hand, the observation of the time delay of the photons in GRBs could provide 
some astrophysical tests for the confirmation of the existence of the extra-dimensions of the 
space-time. 

We have analyzed in detail the time delay for models with both compact and non- 
compact extra-dimensions, deriving some explicit expressions for At. In the case of Kaluza- 
Klein type theories the compactification of extra dimensions provides a natural framework 
for the variation of the fine structure constant and for the speed of light. The basic energy 
scale for this model is of the order of 10 19 GeV, which for Einstein- Yang-Mills type models 
gives a delay that can be as high as At = 10 2 s for 1 TeV photons emitted at a redshift 
of z — 5. In these models the variation in the fine structure constant is dominated by the 
effects of the extra-dimensions, and the quantum effects can be neglected. 

Since in models with large extra dimensions the energy scale can be reduced significantly, 
as a function of the size of v, much higher time delays are expected, which could be of the 
order of 10 5 s for 1 TeV photons. There is a strong model-dependence of the Kaluza-Klein 
type time delay expressions, the results depending on the initial (multi-dimensional) field 
structure. For a pure Einstein gravity in higher dimensions the photon time delay has 
extremely small values, which make it extremely difficult to detect. 

For the detection of quantum gravity or extra-dimensional effects gamma ray bursts offer 
the most reliable high energy photon sources, located at cosmological distances. The BATSE 
data have already been used to extract valuable information on the quantum gravity energy 
scale (Ellis et al. 2003b). However, these data are restricted to a low energy range, of the order 
of 30 — 300 keV. In this range other concurrent physical processes, like Compton scattering, 
could reduce the effect. In order to obtain more reliable data a significant increase in the 
detected photon energy is necessary. The detection of quantum gravitational and extra- 
dimensional features would also require the correlation of GRB redshifts with the temporal 
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and energetic signatures. 

The discovery of the linear polarization of the 7-rays from the GRBs, with the estimated 
degree of polarization of 80 ± 20%, very close to the absolute maximum of 100%, provides 
an other test of quantum gravity effects (Mitrofanov 2003). If the effects of quantum gravity 
are linearly proportional to the ratio E/Eqc, then the polarization of photons with energies 
of about 0.1 MeV should be completely random, contrary to the observations. Consequently, 
quantum gravity effects act with a power greater than one. The linear polarization of 7 rays 
also allow to test the birefringence property of the quantum vacuum, as suggested by the 
quantum gravity in loop representation (Gambini and Pullin 1999). Due to this property, 
two photons with opposite states of helicity have different group velocities (Gleiser and 
Kozameh 2001; Mitrofanov 2003). A significant rotation of the plane of polarization of a 
linearly polarized photon must occur long before any difference in time of arrival is even 
measurable. If a rotation, if present, is below a certain bound, one can obtain a general 
bound on the model parameters characterizing the effect (Gleiser and Kozameh 2001). By 
analyzing the presence of linear polarization in the optical and ultraviolet spectrum of some 
distant sources, the limit x < 5 x 10~ 5 has been found for the dimensionless parameter x 
that characterizes both parity non- conservation and violation of Lorentz invariance (Gleiser 
and Kozameh 2001). This upper bound on x induces a time delay of the order of 10~ 9 s, 
which is beyond the possibility of observation. However, for cosmological GRBs, located 
at a distance of around 10 10 light years, the quadratic birefringence of quantum space-time 
could be tested by polarization measurements of photons with energies greater than 100 
MeV (Mitrofanov 2003). On the other hand, in models with extra-dimensions, there is a 
third order dependence on energy of the speed of light, which could also be tested by using 
polarimetry of 7-rays from cosmological sources. The combinations of both approaches 
based on time-delay measurements and polarimetry could provide significant constraints on 
quantum gravity and multi-dimensional models. 

Ground-based atmospheric Cherenkov telescopes offer a unique opportunity for the ob- 
servation of the delayed TeV components of gamma-ray bursts. In the past few years such 
telescopes using the imaging technique have proved to be remarkably sensitive for the de- 
tection of sources with hard gamma ray spectra (Weekes 2003). Because of the very large 
collection areas associated with these telescopes (> 50,000 m 2 ), they are particularly sen- 
sitive to the detection of transients, e. g. from AGN such as Markarian 421 (Gaidos et al. 
1996). These telescopes are more sensitive than the all-sky viewing ground- level particle 
detectors such as Milagro, the Tibet array or Argo, but have limited fields of view. However, 
they have been used in attempts to detect TeV emission from classical gamma ray bursts 
(Connaughton et al. 1997) and from primordial black holes (Connaughton et al. 1998). Al- 
though the detectable fluence can be as little as 10~ 8 ergs/sec, no detections have been 
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reported yet. However, in no instance has a gamma-ray burst been reported within the field 
of view of an operating imaging atmospheric Cherenkov telescope. Reported observations 
have been limited by the slew time of the telescope, the uncertainty in the initial source 
position and the limited time of operation. 

Therefore, the detection of the time delay between Tev and GeV/MeV photons from 
GRBs could represent a new possibility for the study and understanding of some fundamental 
aspects of the physical laws governing our universe. 
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Fig. 1. — Variation of the photon time delay as a function of the redshift z (in a logarithmic 
scale) in the compactified Kaluza-Klein model, with Einstein- Yang-Mills type action (J3 = 1), 
for 7=1 (corresponding to a fundamental energy scale Ek = 1.2 x 10 19 GeV), and for 
different photon energy values: E 1 = 1 TeV, E 2 = 1 eV (solid curve), E 1 = 10 GeV, E 2 = 1 
MeV (dotted curve), E x = 1 MeV, £ 2 = 1 eV (short dashed curve) and E 1 = 300 keV, 
E 2 = 30 keV (long dashed curve). For the mass and dark energy parameters we have used 
the values Qm = 0.3 and Qa = 0.7, respectively. 
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Fig. 2. — Variation of the photon time delay as a function of the redshift z (in a logarithmic 
scale) in the compactified Kaluza-Klein model with pure Einstein type action (f3 = 3), for 
7 = 10 10 (corresponding to a fundamental energy scale E K = 1.2 x 10 9 GeV), and for different 
photon energy values: E 1 = 1 TeV, E 2 = 1 eV (solid curve), E x = 10 GeV, E 2 = 1 MeV 
(dotted curve), ^ = 1 MeV, E 2 = 1 eV (short dashed curve) and E x = 300 keV, £ 2 = 30 
keV (long dashed curve). For the mass and dark energy parameters we have used the values 
VLm = 0.3 and Q\ = 0.7, respectively. 
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Fig. 3. — Variation of the photon time delay as a function of the redshift z (in a logarithmic 
scale) in cosmological models with large extra-dimensions for a fundamental energy scale 
E F = 7 x 10 15 GeV and for different photon energy values: E\ — 1 TeV, E 2 — 1 eV (solid 
curve), £?i = 10 GeV, E 2 = 1 MeV (dotted curve), E 1 = 1 MeV, E 2 = 1 eV (short dashed 
curve) and E 1 = 300 keV, -E 2 = 30 keV (long dashed curve). For the mass, dark energy and 
dark radiation parameters we have used the values Qm = 0.3, Qa = 0.68 and Qu = 0.02, 
respectively. 
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Fig. 4. — Comparison, in arbitrary units, of the initial KeV/MeV time profile of the GRB 
emission occurring at a redshift z = 3 (assumed to have a Gaussian form), with a duration of 
t = 1 s (solid curve), and the TeV time profile at the detector, modified due to the presence 
of mult i- dimensional and/or quantum gravity effects (dashed curve). Both time profiles have 
been normalized to 1. For the photon time delay we have adopted the linear model, with 
the fundamental energy scale E K = 1.2 x 10 19 GeV. The photon energies are in the range 
Eg (1, 10) TeV. For the mass and dark energy parameters we have used the values Vt M = 0.3 
and fl\ = 0.7, respectively. 



